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ABSTRACT 

We present the angular autocorrelation function of 2603 Dust-Obscured Galaxies (DOGs) in the Bootes field 
of the NOAO Deep Wide-Field Survey. DOGs are red, obscured galaxies, defined as having R - [24] > 14 
{Fia/Fr > 1000). Spectroscopy indicates that they are located at 1.5 < z < 2.5. We find strong clustering, 
with ro = 7.40^0 14 Mpc for the full F24 > 0.3 mJy sample. The clustering and space density of the DOGs 
are consistent with those of sub-mm galaxies, suggestive of a connection between these populations. We 
find evidence for luminosity-dependent clustering, with the correlation length increasing to ro = 12.97^2 64 
Mpc for brighter (F24 > 0.6 mJy) DOGs. Bright DOGs also reside in richer environments than fainter ones, 
suggesting these subsamples may not be drawn from the same parent population. The clustering amplitudes 
imply average halo masses of logM = ILI'^-^Mq for the full DOG sample, rising to logM = 13.0;!;jj 3M0 for 
brighter DOGs. In a biased structure formation scenario, the full DOG sample will, on average, evolve into 
~ 3L* present-day galaxies, whereas the most luminous DOGs may evolve into brightest cluster galaxies. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: statistics — galaxies: high-redshift 
— large-scale structure of the universe 



1. INTRODUCTION 

The bulk of the stellar mass in the universe is created at 
1 < z < 3 (e.g., Dickinson et al. 2003; Rudnick et al. 2006). 
At z w 1 this enhanced star formation occurs primarily in Lu- 
minous Infrared Galaxies (LIRGs; 10" < LirCL©) < 10'^) 
(Le Floc'h et al. 2005), and by z « 2 LIRGs and Ultralumi- 
nous Infrared Galaxies (ULIRGs; Lir(L0) > 10'^) dominate 
the star formation rate (SFR) budget (e.g., Caputi et al. 2007). 
Studies of the spatial distribution of subsets of the 1 < z < 3 
ULIRG population with red optical to mid-IR colors have 
found very strong clustering (Farrah et al. 2006; Maglioc- 
chetti et al. 2008), spanning the range seen from z 2 sub- 
mm galaxies (SMGs; Blain et al. 2004) to high redshift z > 1 
clusters (Brodwin et al. 2007). 

Dey et al. (2008, hereafter D08; see also Fiore et al. 2008) 
presented a sample of IR-luminous Dust-Obscured Galax- 
ies (DOGs) selected via a simple optic al/mid-IR color cut. 
The space density and redshift distribution of these DOGs 
are similar to those of sub-mm selected galaxies (Chapman 
et al. 2005; Coppin et al. 2006). Studies of the spectral en- 
ergy distributions (SEDs) of DOG samples show they contain 
both starburst- and AGN-dominated galaxies, with the AGN 
fraction increasing with luminosity. Bright {F24 « 1 mJy) 
DOGs in Bootes have SEDs of warm AGN ULIRGs (Tyler 
et al. 2008), whereas the majority of faint {F24 > 0.1 mJy; 
(^24)5 = 0.18 mJy) DOGs in the GOODS-N field are domi- 
nated by star formation (Pope et al. 2008b, hereafter P08). In 
this Letter, we study the clustering and environment of DOGs 
as a function of luminosity to explore their role in galaxy for- 
mation at the key z 2 epoch. 

We use a concordance cosmology with 17^ = 0.3 and Ha = 
0.7. Magnitudes are Vega-relative. We report correlation 
lengths in units of comoving Mpc, with Hq = lOOh km 
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s~' Mpc"'. All other physical quantities assume h = 0.7. 

2. DUST-OBSCURED GALAXIES 

We study the DOG sample presented by D08. DOGs were 
identified in the Bootes field of the NOAO Deep Wide-Field 
Survey (NDWFS^, Jannuzi and Dey 1999) via a simple op- 
tical/infrared color selection: R-[24] > 14, or equivalently, 
F24/FR > 1000. Down to a flux density limit of F24 > 0.3 
mJy (« 6a) 2603 sources satisfy this criterion in 8. 140 deg^ in 
Bootes. Spectroscopy (Brand et al. 2007; Desai et al. 2008a,b; 
Houck et al. 2005; Weedman et al. 2005) of 86 DOGs indicate 
that they lie in a relatively narrow range of redshifts, well- 
parametrized by a Gaussian with z = 1 .99 and cr = 0.45 (D08). 

3. CLUSTERING OF DOGS 

The angular autocorrelation function (ACE) of DOGs is 
computed as a function of apparent brightness. Given the nar- 
row range of redshifts (Figure 7 of D08), this binning by flux 
density is to good approximation a probe of the luminosity 
dependence of the clustering. 

The ACF, parametrized as a simple power law, u!(6) = 
Ai^6~^, can be deprojected (Limber 1954) to yield a measure- 
ment of the real-space correlation length, ro(z), over the red- 
shift range spanned by the 2-D sample: 
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Here ■y = 1+S, 
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r(i/2)r[(7-i)/2]/r(7/2), a^(z) is 
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the redshift distribution, and E{z) and x(z) describe the evo- 
lution of the Hubble parameter and the comoving radial dis- 
tance, respectively. The primary uncertainty in the inferred 
real-space correlation length comes from uncertainty in the 
shape of the redshift distribution. 
We calculate the ACF using the Hamilton (1993) estimator: 

DD^RR . 

luh = 1 (2) 

DRxDR 

^ See also http://www.noao.edu/noaodeep/ 
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Table 1. Clustering of DOGs in Bootes 



F24 (^24) 






Halo Mass 


(mjy) (mjy) N 


{hr^ Mpc) 


b 


log(M/M0) L/L,(z = 0) 



>0.3 
0.3-0.5 
>0.4 
>0.5 
>0.6 



0.40 

0.36 
0.53 
0.65 
0.85 



2603 

1846 
1285 
757 
454 



7.40±' 



1.27 
'-0.84 

7.99±lji 

10.19±|^^ 
12.97±«5 



3.12±" 



i_0.51 
■=0.34 
3.36±»| 

3.63±ji;| 



^^■^=■=0.2 
12.3±!!? 



12.5±! 



^^■'=■=0.5 

13.0±0j 



3 4+0'5 
3.7±[|| 

5.0±ii 



Note. — F24 is the 24/im flux density, N refers to the total 
numbers of sources in each flux bin, and b is the hnear bias. Masses 
and luminosities assume h = 0.7. 

^Correlation length from the Bayesian fit. The uncertainty range corre- 
sponds to the 68% confidence interval. 



Here DD, DR and RR are the sum of ordered data-data, data- 
random, and random-random pairs at each angular separa- 
tion. We used 500,000 randoms to ensure a robust Monte 
Carlo integration. We also computed the ACF using the 
Landy and Szalay (1993) estimator and find nearly identical 
results. 

Regions of the survey affected by cosmetic artifacts or data 
quahty issues can compromise a robust measurement of clus- 
tering. Masking is implemented in the images and random 
catalogs to reject these areas. 

A Bayesian technique is used to determine the correlation 
lengths. This allows marginalization over the slope, 6, sub- 
ject to the weak prior that 0.2 < 5 < 1.8. This is desirable 
since the small size of some of the subsamples precludes pre- 
cise simultaneous measurements of both the amplitude and 
the slope. We show in Figure 1 the simple fits and the 
Bayesian likelihood functions in ro computed using the red- 
shift distribution from DOS. The slopes in the fits for all 
samples are consistent with 6 = 0.9. Results are summarized 
in Table 1. 

The ACF errors were derived using the full covariance ma- 
trix, computed using the Brown et al. (2008) implementation 
of the Eisenstein and Zaldarriaga (2001) analytic approxima- 
tion. Conservatively assuming 5% of the DOG sample is spu- 
rious and uncorrelated, then at most we are underestimating 
the clustering by « 11%. 

We show in Figure 2 the correlation length vs. median flux 
density for several DOG subsamples. For comparison we also 
show the 1 (7 range for the Blain et al. (2004) SMG sam- 
ple. The trend is strongly suggestive of luminosity-dependent 
DOG clustering which, as a consequence, implies that bright 
and faint DOGs must reside in different environments. We 
show below (§4) that this is observed. While the uncertain- 
ties preclude ruling out the null hypothesis of no luminosity 
dependence from the clustering measurements alone, the clus- 
tering and environmental studies, taken together, provide con- 
vincing evidence of an intrinsic difference between bright and 
faint DOGs. 

The primary uncertainty in Limber inversion is the sample 
redshift distribution. If the assumed distribution is broader 
than the true distribution the resulting amplitudes will be bi- 
ased high. Our DOG redshift distribution is based on 86 
spectroscopic redshifts, drawn largely from the bright end of 
the DOG population, and may not be representative of the 
fainter DOGs. If the fainter DOGs are largely star-formation 
dominated (e.g., D08; P08) they would likely have a nar- 
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Fig. 1. — Angular correlation functions for (a) the full DOG sample, (b) 
fainter DOGs with 0.3 < F24 (mJy) < 0.5, and brighter DOGs with (c) F24 > 
0.4, (d) ^24 > 0.5, and (e) F24 > 0.6 mJy. The red lines show the best 
fits; the insets show the Bayesian likehhood functions in ro. All slopes are 
consistent with S = 0.9. 
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Fig. 2. — Dependence of correlation length on median flux for the samples 
shown in Figure 1 and summarized in Table 1, illustrating that clustering 
strength increases with median luminosity. The shaded region shows the 1 a 
range for SMGs (Blain et al. 2004). 



tering of brighter DOGs implies they are located in rare, rich 
environments. Indeed, Galametz et al. (2008) find a strong in- 
crease in the incidence of AGN in rich galaxy clusters at z > 1 . 
We show in Figure 3 the surface density profiles of 4.5 fim se- 
lected galaxies from the IRAC Shallow Survey (Brodwin et al. 
2006; Eisenhardt et al. 2004) around DOG samples with sev- 
eral flux limits. We only consider IRAC galaxies with colors 
redder than [3.6] - [4.5] > 0.6, a criterion that selects objects 
at z > 1 .5 (Papovich 2008; Stem et al. 2005). Following Pad- 
manabhan et al. (2008) the mean space density of these 4.5 fim 
sources has been subtracted. All DOG samples are clearly 
correlated with the red IRAC galaxies, showing a large excess 
on small scales. The mean surface densities of IRAC sources 
in the vicinity of the DOGs increases monotonically with their 
brightnesses, indicating that brighter DOGs do in fact reside 
in richer environments, as suggested by their stronger clus- 
tering. The environment of the DOGs is particularly rich on 
small (< 250 kpc) scales, suggesting that they preferentially 
reside in groups of galaxies. This inference is supported by 
their large clustering amplitudes, as well as by recent theo- 
retical work (e.g. Hopkins et al. 2008) showing that at z ~ 2 
the maximal merging efficiency of gas-rich halos, and hence 
resultant starburst activity, occurs in group-mass halos. 



rower redshift distribution due to the strong 7.7/im PAH emis- 
sion feature passing through the 24um filter at z = 2. In 
this case the evidence for luminosity-dependent clustering 
would be strengthened since the correct correlation length for 
fainter DOGs would be reduced relative to that of the brighter 
ones. Quantitatively, if the true redshift distribution were 10% 
(30%) narrower the correlation lengths would be reduced by 
5% (17%). 

4. DISCUSSION & CONCLUSIONS 

To similar flux limits (F24 > 0.4 mJy) the clustering ampli- 
tudes measured for other z ~ 2 ULIRG samples are higher 
than those for DOGs in Bootes. In the SWIRE (Lonsdale 
et al. 2003) survey fields Farrah et al. (2006) find corre- 
lation lengths of ro = 9.4 ± 2.24/1"' Mpc and ro = 14.4 ± 
1 .99/!"' Mpc for ULIRG samples at 1 .5 < z < 2 and 2 < z < 3, 
respectively. In a 0.7 deg- subset of the same survey Maglioc- 
chetti et al. (2008) find r,, = 15.9!^;^/!"' Mpc for a sample 
of 210 z ~ 2 ULIRGs. Beyond the fact that the samples in 
these works have different selection criteria, from each other 
and from the DOG sample, the most likely explanation for 
the differences with our work is that they adopt broader red- 
shift distributions based on photometric redshifts. Desai et al. 
(2008a) and D08 demonstrate that optical/IR photometric red- 
shifts can be unreliable for these heavily obscured sources, 
particularly when the optical detections are marginal or non- 
existent. These previous analyses have likely underestimated 
their uncertainties, because they fixed the slope of the corre- 
lation function, eliminating the covariance between the slope 
and ro, and adopted simple Poisson errors, ignoring correla- 
tions between adjacent bins. 

While DOGs are a mixed population, consisting of both 
starbursting galaxies and AGN, a majority are dominated by 
star formation (P08). It is interesting then that the space den- 
sities (D08) and clustering are quite similar to SMGs (Blain 
et al. 2004; Coppin et al. 2006), which are known to be star- 
formation dominated (Pope et al. 2008a). 

On the other hand, the AGN fraction increases with lumi- 
nosity (P08; D08; Tyler et al. 2008). The very strong clus- 



'E 

3. 

03 
X 

a 

CB 

(5 



E 0.6 



>< 0.4 



0.2 

03 
£= 
03 

Q 



03 

O 

CO 

in 
w 

03 
O 
X 
LU 



0.0 



-0.2 



Radius (Mpc) 
0.0 0.5 1.0 1.5 2.0 2.5 



• F24 (mJy) > 0.6 
o F24 (mJy) > 0.5 
A 0.3 < F24 (mJy) < 0.4 



^ /| A jil # ih 



1 2 3 4 5 
Radius (arcmin) 



Fig. 3. — Surface profile of z > 1.5 4.5/.tm sources in Bootes around DOG 
samples with various flux limits. The contribution from the mean has been 
subtracted. There is a clear luminosity dependence to the profiles such that 
brighter DOGs lie in increasingly rich environments. The top axis is in phys- 
ical separations at z = 2 for our chosen cosmology. The region within < 250 
kpc of the DOGs is pailicularly rich. Some of the symbols are offset slightly 
in radius for clarity. 



D08 proposed that both types of DOGs are drawn from 
the same parent population, and that the luminosity depen- 
dence of the AGN fraction arises from the additional energy 
output from those DOGs undergoing an active AGN phase. 
The present results suggest an alternative explanation. Al- 
though the evidence for luminosity-dependent clustering is 
marginal given the large errors, the corroborating observation 
that brighter DOGs reside in richer environments than fainter 
ones indicates that they are not drawn from identical parent 
populations. This conclusion is robust to uncertainties in the 
redshift distribution, provided a single distribution is used for 
both bright and faint DOG samples. If the redshift distribu- 
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tion of faint DOGs were narrower, the likeliest situation, the 
strength of the luminosity dependence would increase. 

The Unear biases of the DOG samples, Usted in Table 1, 
are computed as the square root of the ratio of the z = 2 
DOG and dark matter correlation functions at a scale of 5 
Mpc, where the latter is computed following the HaloFit 
prescription of Smith et al. (2003). The full sample has a 
bias ofb = 3.l2'^'ll, where the uncertainty is propagated from 
the error in the clustering. The bias increases with DOG 
flux, from b = 3.36!*]:| for faint DOGs ((F24) = 0.36 mJy) 
tob = 5.33^1:^1 for bright DOGs ((F24) = 0.85 mJy). Compar- 
ison of the observed clustering with that of halos in a large, 
high-resolution numerical simulation (described in detail in 
Brown et al. 2008) indicates that the full DOG sample has 
an average halo mass of logM = II.T^^'IMq. The masses 
increase with DOG luminosity, as shown in Table 1, reach- 
ing logM = 13.0!jJjMQ for the brightest DOGs. In the Fry 
(1996) biased structure formation model, assuming merger- 
free passive evolution, these samples will evolve into « 3.4 
and sa 6.6L* galaxies by the present day. The latter have the 
masses of brightest cluster galaxies in local clusters. 

At low redshift the space density and color bimodality of 
galaxies can be modeled by truncating star formation in galax- 
ies above a particular host halo mass (e.g., Brown et al. 2008; 
Croton et al. 2006; Kauffmann et al. 2003). Several papers 
suggest that the transition halo mass is ~ lO'^ M©, and this 
mass undergoes negligible evolution at z < 1 (e.g.. Brown 
et al. 2008; Cattaneo et al. 2008; Dekel and Birnboim 2006). 
The z ~ 2 DOGs, most of which are undergoing vigorous star 
formation, reside in w 10^^- 10^^ M0 halos. Perhaps the 
mode of gas accretion onto massive halos changes at z > 1, 
as suggested by the virial shock heating model of Dekel and 
Birnboim (2006). 
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forming galaxies with halo masses larger than lO'^M©, the 
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They will likely evolve into very massive (3 < L/L* < 7) lo- 
cal galaxies. 
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